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ABSTRACT 


Orbiting astronomical observations have the potential for making 
observations far superior to those from earth -based mirrors. In order 
for this performance to be realized, the contour of the primary mirror 
must be very accurately controlled. A thermally activated system for 
correcting symmetrical distortions in space telescope mirrors has been 
evaluated. This system utilizes thermally induced elastic strains to 
correct axial distortions in the mirror. The relation between axial 
distortion and thermal inputs was determined hy a finite difference 
solution of the equations for thin elastic shells. 

The use of this technique was demonstrated analytically on a 
beryllium paraboloid. This mirror had 10 equally spaced thermal inputs. 
Distortions due to an acceleration-type loading were shown to be corrected 
to well within the required accuracy. Axial temperature gradients 
resulting from the application of the thermal inputs to the rear surface 
of the mirror were shown to be quite small. 
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CHAPTER I 


INTRODUCTION 

The resolution of ground-hased optical astronomical observa- 
tories is limited by atmospheric turbulence. To minimize the effects 
of turbulence, many observatories are located at relatively inaccess- 
ible areas atop tall mountains. Even at these locations the very best 
telescopes seldom have a resolution better than 0-3 second of arc. 

In space, telescopes would he limited only by the diffraction limit 
and, therefore, large mirrors (l$0 in. dia. ) should he able to resolve 
0.03 arc second at 5>000 angstroms (Ref. l). Also, space telescopes 
will be able to view portions of the spectrum not currently available 
for observation due to atmospheric absorption. 

The National Aeronautics and Space Administration is currently 
investigating the problems associated with the operation of large 
(120 in. dia. aperture) space telescopes. Several studies of space 
telescopes have been conducted and are reported in References 2, 3> 
and k- These studies have defined the scientific objectives of a 
space telescope and outlined some of the major problems involved in 
its design and fabrication. All of these studies have recommended 
the use of a Cassegrain optical system which requires a paraboloidal 
primary mirror. In order for the optical system to have the desired 
resolution, the contour of the primary mirror must he maintained to 
within 2 microinches of the design value (Ref. 3)* If this accuracy 
cannot be maintained, the resolution of space telescopes will be 
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degraded relative to their potential capabilities and their performance 
may be less than that of ground -based observatories. 

Two types of primary mirrors have been proposed for use in space 
telescopes. The first type is a passive mirror which would be designed 
to retain the proper contour without correction for the life of the 
telescope. This mirror would consequently he rather massive and may 
impose a severe weight penalty on the launch booster. The second con- 
cept is known as an active optics system and utilizes a thin mirror 
which is peimitted to deform moderately under operational loads. 
Distortions in the mirror would be monitored and analyzed by a figure 
error sensor. This sensor would activate a control system to apply 
correction loads which would remove the distortions. 

One active optics system utilizing precision jacks to provide 
corrective loads has been investigated analytically. This analysis has 
been experimentally verified using a thin deformable mirror 30 inches 
in diameter (Ref. 5). However, this system requires either a very 
stiff back plate for the jacks to react against or a determination of 
the coupling between the hack plate and the mirror. Also, this is an 
electromechanical system and is relatively complicated for space use. 
Another type of active optics system that has been suggested utilizes 
thermal inputs to provide the corrective distortions. Elastic strains 
introduced hy differential heating would be used to force the mirror 
to assume the proper contour. 


The object of this investigation is to develop a technique for 
determining the relation between deformations parallel to the mirror 
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axis and the thermal inputs necessary for error correction in such a 
system- Thermal inputs in the form of a prescribed temperature distri- 
bution were considered to be applied to the rear surface of the mirror. 
In order tq illustrate the feasibility of such a technique, a cosine- 
type temperature profile was considered for the loading. The axial 
temperature gradient introduced in the mirror due to front surface 
radiation loss wps shown to be small. The relation between the mirror 
distortions and the thermal inputs was obtained by a computerized finite 
difference solution of the elastic shell equations. The relation was 
expressed in the form of a flexibility matrix. The thermal inputs 
necessary to correct distqrtion at specific control points were deter- 
mined by inverting the flexibility matrix to form a stiffness matrix. 

An example of t^e thermal inputs necessary to correct distortions due 
to ap acceleration-type loading is shown. 



CHAPTER 11 


THERMAL CONTROL SYSTEM 
description of Telescope 

The space-telescope model selected for this analysis is a pre- 
liminary design of the type discussed in Reference 3- -A sketch of 

this model is shorn in Figure 1. The "basic telescope configuration 

* 

consists of two large cylindrical shells which are attached to the 
telescope cabin. The cylindrical shells enclose the main optical 
elements - the primary and secondary mirrors. ; Attached to the outer 
shell is a system of doors that prevents sunlight from falling on the 
optical system during maneuvers. The 'inner shell is a thermal, shield 
which reduces solar heating loads on the primary mirror. All optical 
imaging devices and sensing instruments are contained in the telescope 
cabin. This cabin will 'also provide the necessary environment for ’ 

4 

manned support. One significant departure from previously designed 
space systems is that the telescope must be capable of continuous 
operation for several years. Therefore* the telescope must receive 
manned support from a docked or nearby space station. 

The telescope has a Cassegrain -type optical system with a 
focal -length-to-aperfcure diameter (f/d) of 30* The primary mirror is 
a short -focal-length (f/d = 4 ) paraboloid, while the smaller secondary 
mirror is a, hyperboloid. If this optical system could he fabricated 
with perfect geometry it would image a point source, such as a star, 
in. the focal plane as a bright central disc (Airy disc) with 
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surrounding diffraction rings. The Airy disc would contain approxi- 
mately 85 per cent of the incident energy with the remaining 15 per cent 
scattered in the surrounding rings. However, due to fabrication errors 
in the mirrors, operating loads and inherent absorption losses, the 
energy in the central disc will be reduced and the resolving power of 
the telescope may be severely limited. If the distortions in the 
optical system are such that 68 per cent of the energy lies in the 
central disc with 32 per cent in' the diffraction rings, the 'system is 

k 1 

\ , ! 

operating at the Rayleigh criterion for resolution. An analysis of the 

1 . * 
optical system for this diffraction-limited operation was cbnsidered in 

* 1 ' * - k 

* 

Reference 3. This analysis estimates that, the primary mirror will • < 

■»- ‘ , ■: • ■ 

require a surface contour having distortions less than 1.97 microinches 
(l/lO wavelength at 5*000 angstroms) from .the design paraboloid. In 
addition, a ro ot -me an - s quare surface accuracy of 0-37 microinch (1/53 
wavelength at 5*000 angstroms) will he required. 

'H _ i 

It was noted previously that distortions in the optical system 
may he produced by operational loads. Other sources of distortion may 
include (l) the introduction of elastic strains in changing from an 
earth gravity environment to a zero gravity space environment, 

(2) relaxation of residual strains introduced during the fabrication of 
the mirror billet, (3) relaxation of residual strains- introduced by 
machining and polishing operations, and ( 4 ) plastic strains introduced 
by launch and environmental loads. 



Thermal Active Optics System 

In order to maintain ' the mirror contour within the required 


accuracy , a technique to introduce corrective ' distortions may he 
necessary* The technique envisioned in this analysis for correcting 
mirror distortions is to apply thermal inputs to the rear surface of 
the mirror* These thermal loads will induqe elastic strains which will 
deform the mirror surface to the desired, paraboloidal contour. A 
schematic diagram of a thermal active optics system is shown in 
figure 2. Distortions in the primary mirrpr are detected by the figure 
error sensor using interferometric techniques. These errors are in the 
form of fringe patterns apd must he interpreted to determine the size 
^nd direction of the distortion. This interpretation is performed by 
the analyzer and phase detector at fixed control points an the mirror 
surface, The analyzer will also calculate the amplitude of the con~ 
trolled temperature source necessary to correct distortions at the 
control points* The function of the control system is to apply the 
desired inputs, thereby reducing the distortions at the control points 
to (or below) the acceptable level. 

A sketch illustrating the control points and location of the 
thermal inputs is shown in Figure 3 * The thermal inputs are in the 
form of controlled back surface temperature distributions. The radial 
and circumferential location of the thermal inputs, in addition to 
their size, may vary. It is only necessary to determine the Influence 
pf the thermal inputs at specific control points on the surface. 
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Figure 3- - Sketch showing control points and location of 
thermal inputs on rear surface of telescope mirror. 

Most practical structures, including the telescope mirror dis- 
cussed in Reference 5 have teen shown to exhibit approximate linear 
behavior. In order to determine the influence of the thermal inputs, 
it was assumed in this analysis that the thin paraboloidal shell 
behaves as a linear structure. Thus the influence of thermal inputs 
can be analyzed using linear theory. One advantage of such a struc- 
tural theory is that it permits the application of the principle of 
superposition. This principle states that stresses and deformations 
produced in a structure by a set of loads in combination can be obtained 
by adding the stresses and deformations produced hy each load acting 



10 

separately. Therefore, the distortions at the control points may he 
expressed as a function of the applied thermal inputs by the following 
equations : 

6 1 = c ir!rl * C 12 T 2 + c 13 T 3 + c l4 T 4- + * * * c ln T n 

S 2 = c 21~l + c 22 T 2 + c 23 T 3 + c 2b\ + * * * c 2n T n 

* 
m 

5 n ~ c nl~l “ G n2^2 "* Q n3 T 3 + 0 'n4 T 4 * * * * c nn^n 

t , 

The coefficients ( ) specify the contribution due to an increase in 
temperature from, each thermal input (T^ ) toward the distortion (S^ ) at 
the control point (i). These equations may be written in the more con- 
venient matrix notation as 

1*1 - [<G hi (o 

The square matrix fcj is generally called the structure flexibility 
matrix and the component terms are deflection influence coefficients. 

The thermal inputs necessary to correct mirror distortions can 
be expressed by inverting the flexibility matrix and multiplying by the 
measured distortions. 

M-BM 0) 

DG - DO " 1 



where 
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In order for the active optics system to perform, properly it is 
necessary to accurately determine the coefficients k.^ for each 
control point and store the coefficients in the system analyzer. 



CHAPTER III 


THERMAL ANALYSIS 
Description of Mirror Geometry- 

No firm design of the space telescope has "been formulae ed. The 
technique developed in this investigation is applicable to any thin 
telescope mirror. The example shown in this and subsequent sections 
illustrates the use of this technique. 

A sketch of the mirror used in this analysis is shown in 
Figure 4-. The mirror is a thin paraboloidal shell having a diameter 
of 120 inches and a focal length of 480 inches. Since no firm design 
of the telescope has been formulated, certain assumptions concerning 
the mirror geometry were necessary. For example, the diameter of the 
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center hole mil depend upon the telescope optical system and manu- 
facturing considerations. The hole was assumed to have a diameter equal 
to one-tenth (12 inches) the diameter of the mirror. In the weightless 
environment of space, only a very thin reflective surface would he 
necessary for the primary mirror. However, practical considerations 
of manufacture will require the mirror to he sufficiently thick to with- 
stand grinding and polishing in a gravity environment. For this invest! 
gation, the thickness of the mirror was one-hundredth (1.2 inches) of 
the diameter (Ref. 7)* The thickness was assumed to he constant in 
both the circumferential and meridional directions. 

Metal mirrors are ideally suibed for thin one-piece construction 
because they have high stif f ne s s -t o-w eight ratios. Beryllium has one of 
the highest stiffness -to -weight ratios of any structural metal nnd is 
currently being considered as one of the prime candidate mate <’j.als for 
telescope mirrors (Ref . 8). The paraboloid shown in Figure *i was con- 
sidered to he fabricated mom a homoge.jenu.i and isotropic billet of 
beryllium. The properties of the beryllium material for this mirror 
are shown in Table I. 

The mirror was restrained at the outer rim by a system which 

accommodates only symmetrical loading about the mirror axis. This 

support system restrains the mirror only in the axial direction and is 

usu ally referred to as a hinged support on rollers. A sketch of the 

* , ' ’ 
support condition is shorn ip Filgure 5 . , This system is similar to a 

three-point tangent-bar mounting suspension considered in Reference 3 . 

The systems are similar in that both will accommodate differential 
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Figure 5° “ Sketch of mirpor cross section illustrating 
boundary conditions. 

radial expansion between the mirror and support structure without 
introducing loads in the mirror- Ifo restraints were’ applied at the 
central cutout portion of the mirror. 

For the illustration of the thermal active optics system, 
symmetrical distortions of the primary mirror were considered. Sources 
with a controlled temperature distribution were selected to provide the 
thermal inputs. The controlled sources were applied to the rear 
surface of the mirror. In ord,er to simplify the problem, only steady- 
state distortions were examined. Therefore, the results of this study 
are applicable only if observations are made after steady-state 


conditions exist. 



Ten equally spaced stations along the meridian of the mirror 
were selected as the control points. The location of these stations 
is shown in Figure 6. The points wqre separated by a meridional dis- 
tance of 5 .!jC 4 inches. Ten control points were considered to be a 

sufficient number to demonstrate this technique. The use of more 

1 * 

control points would only have generated a larger matrix of influence 
coefficients and would have added little to the demonstration of the 

S 

technique. For an actual control system, a larger number of heaters 

may be desirable to increase the control capability. No cpntrol point 

* ' 

was located at the periphery of the shej.1 because of the axial restraint 
imposed by the boundary condition. A thermal input at that location, 
however, would influence the distortion at the other control points. 

The thermal inputs applied by the strip heaters to the rear 
surface of the mirror were assumed to have the form 


M 


; + 1 cos 

V 2 ' 2 t ' 




(5) 


= 0 
v 


I* - 


> 


A 


A sketch of this thermal input applied at a control point on the mirror 
surface is shown, in Figure J. Thermal inputs of this form were selected 
because they can be represented by a concise mathematical expression. 
Also, these inputs should be relatively easy to simulate experi- 
mentally since they produce no severe radial gradients. It should 
be noted in the figure that the thermal inputs were applied to one-half 






18 


of the interval represented by the control point. This permits a 
spacing between the strip heater?. 

The thermal inputs of Figure 7 were used at every control point 
except station 1 (see Fig. 6 ). At station 1 a slight modification of 
the thermal ipput was psed. Only the right half of the distribution 
shown in Figure 7 was applied and the interior of the hole was assumed 
to be insulated. 


Internal Temperature Distribution 
The temperature distribution indicated ip equation ( 5 ) was 
applied to the rear surface of the mirror at each control point. The 
application of this axially symmetric thermal input will result in 
two-dimensional heat flow within the mirror interior. Therefore, the 
interior temperature of the mirror (T^) will be 

( 6 > 

where £ is the meridional coordinate and £ is the normal to the 
neutral surface. This temperature distribution must be determined in 
order to evaluate the effect of the thermal inputs in reducing mirror 
distortions. The internal temperature distribution could be signifi- 
cantly affected by the radiation heat loss from the mirror front sur- 
face. In order to evaluate the front surface heat loss at each control 
station, a thermal model of the telescope was examined. This model is 
shown in Figure 8 . The mirror was' represented hy a flat circular disc. 

A disc was considered to be a good approximation because the mirror is 




VO 


Figure 8.- Thermal model for calculation of axial heat loss. 
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a very shallow paraboloid with an. edge to center depth of approxi- 
mately 1.9 inches. The temperature of the mirror and thermal shield 
are 407° R and 390° R* respectively, and are assumed to be constant 
over the surface. These temperatures are based on a preliminary heat- 
transfer analysis reported in Reference 3 . The interior of the thermal 
shield was coated with lampblack which has an emissivity of 0.96. The 
purpose of this coating is to eliminate stray radiation from falling 
on the telescppe optics and introducing noise in the observations. 

Also, the secondary mirror and support struts were omitted from the 
model because the area was considered small. 

The front surface radiation emitted by the mirror will he 

~ c ^^V^m (?) 


where the emissivity and absorptivity are considered equal since the 
mirror and thermal shield are at approximately the same temperature. 
The mirror temperature is not 'considered to be significantly increased 
by the thermal inputs. The radiation from the thermal shield that 
strikes the mirror is 


Q s-*T a s a 


T F 
s s sm 


( 8 ) 


where 3?^ is the view factor and represents that portion of the 

energy emitted by the thermal shield which is intercepted by the 

1 

mirror. Neglecting reflections, the pet heat loss from the front 
surface of the mirror is equal to the radiation emitted by the mirror 
( oa^A^T^ ) minus the radiation absorbed by the mirror from the thermal 
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shield (aa s T s ^A s F sm • a^). This is expressed in equation form by the 
following relation: 

Q — crajjjA^Tjjj - aa s T s •^■s^’sjn a m ( 9 ) 

In or^er to calculate the heat loss from the mirror wi^l be necessary 
to evaluate the view factor (F sm ) relating the radiant emission from the 
thermal shield which goes directly to the mirror. Since view factors 
are generally difficult to calculate due to involved integrals, simpli- 
fying relations are often sought. In order to determine the view faptor 
between the thermal shield and mirror (F gm ), it will he related to the 
view factor between the mirror and open port (F m p) which can he readily 
evaluated. 

Since the mirror, heat shield, and opening at the end of the 
telescope (port) effectively form an enclosure and the disc does not 
radiate to itself. 


F + F = 1 


ms 


mp 


Using this equation and the reciprocity relation 


( 10 ) 


^Wuns (H ) 

the view factor F^ can be determined. By equation (11) 

, , Vms 

Sm A s 


( 12 ) 
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Substituting for from equation (10) 


' W 


Replacing F^ in equation (9) by equation (lj5) yields 

Q “ An cra m|%i^ " ^s^ a s(l ” Fmp)j ( llj 

In order to solve equation (l4 ), it is necessary to determine the view 
factor F m p. A sketch of the geometry is shown in Figure 9* From this 
sketch and the definition of the view factqr (Ref. 9), we have 


■V^mp f f 

^ J *y 


cos COS f p dAp di^ 


where 


e°s = cos f p = - 


h - N/L 2 + (p^ - P p ) 2 (2 

Substituting the above equations into equation (l^ ) and also substitu- 
ting for dAp gives the following equation: 


p p 2 jt p 8o 

vw=j. L L 


^ ° ° *(l 2 + (Pp - p m) 2 “l 


2 ap p d0p 





Integrating this equation over dp^ and &0p and substituting 
L = 430 in. yields the following integral over the mirror area: 



1.84-9 X 10 5 

1.84-9 x 10 5 + (80 - p m )2 


Pm (80 ” Pm) 

1.84-9 X 105 + (80 - P m ) 2 


+ 


4-30 



80 - Pm \ 

v 430 / 



(19) 


Only the heat loss from the front surface at the control points 
is of interest. Eaoh control point represents a small portipn of the 
total mirror area. Therefore, since A m is small in relation to Ap, 
equation (19) gives the view factor F m p directly in terms of the 
mirror radius where it is to be evaluated. Therefore, 


F =1 
mp 


1.849 x 10 5 


P m (8o - P m ) 


1.849 x 10 5 + (80 - p m ) 2 1.849 x 10 5 + (80 - p m ) 2 


+ 




80 ~ p m\ 

^30 I 


+ tan 



( 20 ) 


All terms necessary for evaluating the heat loss (eq. (l4)) are 
now known. The calculated values for the view factor and the heat flux 
(Q,/A m ) at each control point are shown in Table II. The view factors 
are quite small and do not vary significantly with the location of the 
control point. Since the heat flux is directly proportional to the 
view factor, it is also quite small. 
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The temperature gradient at the mirror surface can he determined 
using the Fourier heat conduction equation 

'Ha ^ = Q (21) 


Simplifying this equation, 


dT 1 Q 
dl " K \ 


( 22 ) 


An examination of a typical temperature gradient (control point 5) 
indicates that 


|| = 1.861 X 10 J ' °R/in. (25) 

This is a small temperature gradient due to the view factor and the 
high thermal conductivity of "beryllium. Since the temperature gradient 
at the mirror front surface is small, this boundary may he assumed to he 
insulated (dT/d^ = ^_ - 0) when determining the interior temperature 
distribution. Several nonmetallic materials, including ceramics and 
glasses, are also being considered for mirror fabrication. These 
materials have a very low thermal conductivity and therefore would have 
a larger temperature gradient at the mirror front surface. 

The interior temperature distribution was examined try applying 
a thermal input of unit amplitude to the flat circular disc. A section 
view of the disc illustrating the coordinate system and thermal input 
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is shown in Figure 10. The thermal input is applied over the interval 
P.j. • The front surface is thermally insulated due to the low-temperature 
gradient discussed previously. It should he noted that the conditions 
considered in this analysis are somewhat different from those shown 
in Figure 7 . The thermal input is applied to the circular disc from 
'the center outward, whereas the actual telescope would have the central 
portion removed. The solid disc was chosen because it should have 
little effect on the internal temperature distribution and the finite- 
ness condition at the center (discussed later) readily permits evalua- 
tion of constants necessary for the solution of the differential 
equation. 



Figure 10. - Sketch illustrating the coordinate system and 
boundary conditions for deteimining the temperature 
distribution in a circular disc. 



The general differential equation for an. axisymmetric steady-state 
temperature distribution in a circular disc is given below (Ref. 10): 
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a g T i ar aS 
ap 2 p a? a? 2 


(24) 


The boundary and finiteness conditions applicable to this 
problem are 


T(0,£) is finite 


(25) 


aT 

as 


S— t 


= o 


insulated front surface 


(26) 


Therefore, 


T(p,0) - A| 


1 

2 


1 cos 

2 pj . 


(27) 


ai 

ap p=p t 


= 0 


(28) 


The finiteness condition (eq. (25)) results from physical limitations 
on the temperature at the center of the disc. It was noted that the 
disc was chosen instead of an annular ring because the finiteness 
condition simplifies the solution without significantly modifying the 
problem. The second boundary condition (eq. (26)) has been discussed 
previously. Equation (28) can be shown directly from the imposed 
temperature distribution given as equation (27). The equation for the 
applied distribution (eq. (27)) is different from equation (5 ) i however. 



the actual thermal input is the same. The governing differential 
equation (eq. (24 ) ) can be solved by the separation of variables 
technique. Assume a produqt solution of the fqrm 
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T(p,g) = B(P)2(6) (29) 

Substituting equation (29) into equation (24) and performing the 
indicated differentiation yields the relation 

1 1 ;» _ 19 s ! . . 

Since the left-hand portion of the equation is independent of £ and 

the equivalent right-hand member is independent of p, bpth sides must 

therefore be independent of £ and P and may be set equal to a 
2 

constant -A . 


isV 1 Sr 

"t" — ■ 1 ■ “ ■ " 

R cp Sp 



(3D 



(32) 


Equation (31) is a Bessels equation. The solution to this equation is 
a Bessels function of the second Rind of order zero. Equation (32) is 
an ordinary linear differential equation whose solution may he obtained 
using qperator techniques. The solutions to both equations are given 
in equations (33) and (34), respectively. 
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R = C lt J 0 (Xp) + C 2 Y q (Xp) 


(33) 


Z = C 5 e A ^ + C4 e'^ (34 ) 

The temperature may now be determine! from the product of equa- 
tions ( 33 ) and (34). 


T = R(p)Z(£) 


“ [CiJ 0 (Ap) + a 


,Y 0 (Ap)3 O 7 




C 4 




(35) 


The constants will be determined by imposing Ijbe boundary and finiteness 
conditions. The first condition requires that C 2 = 0 because the 
limit Y 0 (Xp Redefining and combining constants, the temperature 

p 0 

may be expressed as 


T = 




+ C 2 e 




J 0 (XP) 


( 36 ) 


Partial differentiation of equation (36) with respect to £ and 
imposing the boundary condition of equation (£ 6 ) yields the following 
relation between constants G-^ and C 2 where t denotes the mirror 
thickness. 

C 1 = C 2 e ’ 2 ^ (37 ) 


Therefore, 


T = 0 2 J o (Ap) 


~hi"2ht 

e 



(38) 



The separation constant A may be defined by the fourth boundary 
condition given in equation (28). 
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St 

8 p 


p=pt 


= -CgAJ-^Apj. ) 


" At -2At , -AS 
e + e 


= 0 


( 39 ) 


The only way this equation may be satisfied without having a trivial 
solution is for 


AjqCAPt) = 0 


(4o) 


which yields the following first four values for A: 


A q = 0 


A 


P 


3.8317 


Ao = 


t 

7.OI56 


A = 10 ‘ ^3 
5 ~ p t 


j 


(U) 


The general solution of the equation must involve the sum over all A 
and may be written as 


T 


- 2 C 0 + £ C n J 0 (V> 


n=l 


e A n S-2A n t + e “^ n = 


L 


(^2) 


The above equation specifies the interior temperature distribution of a 
flat circular disc with the front surface insulated against heat losses. 
The relation is complete except for determining the constant C which 
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may "be defined using the imposed temperature distribution. By examining 
the distribution over the interval from 0 to p and noting the 
orthogonality relation,, the constants can be determined from equations 
defined in Reference 11. From the reference, the constants are 


C 


o 




to) 


c „ = 


/ Pt pT(o,o)j 0 (v)to 


B n (l + e ) 


B n =/ 1 p C J o ( V0 S 40 


(It) 


to) 


where n = l,2,3,t, . . . 


Substituting a (~ - | cos ~-j for T(p,0) into the equation 


Pt 

for C D , that constant can be found to be 


„ A 

o "1? 


Performing the same substitution, C n may be found to be 


to) 


G n-“ 


2B n (t 

r 

where B n is given in Reference 11 as 


77®^ f 0 1 p ( 00s yf) J ° <7 '“ p)dp m) 




m 
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Therefore, substituting C Q into equation (42) and 3^ into equa- 
tion (47), we have 


00 



The integral in equation (50) for C n cannot be evaluated directly 
but must be found graphically or numerically for each value of 
A computer program to evaluate both equations (49) and (50) was written 
in the Fortran 2.3 programing language and is included in Appendix A. 

The Gauss quatrature method was used to evaluate the integral in 
equation (50)* The computer program was utilized to determine the 
difference between the applied ba,ck surface temperature and the tempera- 
ture of the interior as a function of radius. The first 20 terms in 
the series of equation (49) were used. 

The results of these calculations for a unit amplitude indicate 
that the maximum temperature difference is 0.0l6° and occurs along the 
axis of the disc. The temperature difference as a function of radius 
for £ - 0.5 and £ - 1.0 is shown in Figure 11. It should be noted 
that the maximum temperature difference does not occur at the same 
location as the maximum thermal input. 




Figure 11. - Temperature difference through a circular disc as a 

function of radius. 


CHAPTER IV 


THERMAL 3HPLUEHCE COEFFICIENTS 

i 

Paraboloid of Revolution 

The distortion influence coefficients relating the thermal input 
and the distortion at the control points w ill now be determined. Several 
techniques for determining these coefficients are available. One method 
which has been employed for the force active optics system is the finite 
element technique (Ref. 12). This technique consists of dividing the 
mirror into small interconnected elements of finite size. The deforma- 
tions of the mirror are determined at the points of connection called 
nodes. The method used in the present investigation is the finite 
difference technique in which the governing differential equations are 
solved by approximating the derivatives by finite differences between 
nodes. Either technique may be used to determine the deflection 
influence coefficients. Once the deflection influence coefficients 
(c - - ineq. (l)) have been evaluated, a flexibility matrix for the 
control points can be formulated. The amplitude of the thermal inputs 
necessary to correct mirror distortions can he expressed by inverting 
the flexibility matrix. 

The equations governing the linear behavior of thin shells of 
revolution are well known and may he found in References 15 and lip. 

The basic equations will be shown here to clarify this analysis. The 
mirror was considered to be fabricated from a homogeneous and isotropic 
material and only static symmetrical distortions are of interest. The 



35 


geometry and coordinate system tor a paraboloidal shell are shown in 
Figure 12. Any point in the shell may be located by specifying its 
coordinates (£ ; 9,£ ). The origin and positive direction of the coordinates 
are indicated on the figure. The meridional coordinate is denoted as I, 
the circumferential coordinate as 6, and the normal to the tangent plane 
is indicated as The neutral surface is chosen so that 

J t E d£ = 0 (51) 

where the integration is through the thickness. This permits a variation 
in the modulus due to temperature changes to be considered. The modulus 
of the beryllium used in this analysis was considered to be constant; 
therefore, the neutral surface will be the middle surface of the shell. 

The principal radii of curvature Rg and Rg are written in terms of P 
and | as 



These relationships can be defined using the usual parabolic relation- 
ship between the radial and axial coordinates. 

( 5*0 


p 2 = 
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Utilizing the definition for the length of a differential element 



+ <3y 3 


(55) 


we can obtain 



(56) 


from which the radii of curvature can be shown to be 


R e 



(57) 


R, 



x (Ij-f 2 + p *) 5 / 2 


(58) 


A shell element indicating the positive directions of the membrane 
forces per unit length, transverse forces per unit length, and loading 
per unit area is shown in Figure 13(a), The moments per unit length 
are given in Figure 13(b) and the positive directions of displacement 
and rotation are shown in Figure Ik. The equilibrium equations for any 
isotropic shell of revolution loaded axi symmetrically are shown below 
(Ref. 15 )• 


3pn, 


" • N « 3s + s^-sr ' 9 sr 


) 




= 0 




S /5pMs 5p 

' Me sr. 


) . 3 + a' 

/ V E i K e/ 


+ Pq = 0 , 


(59) 



(a) Stress resultants and loads. 


(b) Moment resultants. 


Figure 13.- 
moments. 


Figure 14.- 


Shell element with stress resultants, stress 
and loads applied in the positive sense. 



Shell element with displacements and rotation 
indicated in the positive sense. 
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The transverse shear forces have "been eliminated by using the moment 
equilibrium equations. 

The displacements and rotation of the midplane surface are 
relabed by 


i bw u | 

$ t = - -- + Js. 

* E| 

• * i 

The membrane strains are expressed in terms of displacements by 


dug w 

€. s + — 

6 dg Rg 


_ u | dp w 

Se = 7 ¥ + i: 


J 


The distortions due to bending are given by 


K t = 


d| 


cr D - vo t 


(60) 


(61) 


Sp 

Ks "T 5 e 


Neglecting the effects of stresses normal to the shell and 
assuming surface normals to the neutral surface remain normal after 
deformation^ the stress-strain relations are given by the following 
equations : 

Of - V0Q 

H + CKg = + aT 


(62) 


( 63 ) 


•f aT 


e 9 + ^ K e 55 


E 
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The stress- strain relations along with the definition of the stress and 
moment resultants 


h - f°t d£ 


(64) 

h 6 => f u e a? 

Mg =jT £a e d£ 



and the relation 

f & di = 0 (51) 

yield the following relations: 
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Equations (59) - (62) and (65) constitute a set of 11 equations relating 
11 independent variables. Ey combining these equations, a set of three 
second-order differential equations in terms of meridional and axial 
displacements (ug and w) and meridional moment resultant (M| ) can be 
obtained. 

The three second-order differential equations are subject to 
restraints applied along the boundary of the shell. The boundary 

V ♦ 

conditions and mirror support system were discussed previously and are 
shown in Figure 5- The periphery of the shell will not support a 

1 

moment Mg due to the hinge. The sum of forces in the radial direction 

equals zero, as do the axial displacements. The positive direction of . 

forces and displacements at the outer -boundary are shorn in Figure 15(a). 

* 

The equations used in describing this boundary condition are as follows: 

Ms = 0 r 

[■ 

w cos cp - U| sin <p = 0 \ (66) 

sin cp + cos cp = 0 

The central cutout portion of the mirror was. unrestrained. 
Therefore, no resultant forces or moments can be accommodated at this 
boundary. The positive direction of forces and displacements at the 
inner boundaiy are shown in Figure 15(b). The equations used to 
describe this boundary condition are as follows: 





b-3 

Mg = 0 
Ng = 0 
% = 0 

As noted previously, diffract! onr-limi ted operation of a telescope mirror 
requires that axial distortions of the primary mirror be kept within a 
tolerance limit of 2 microinches. The relation between axial distor- 
tions and the distortions along the shell meridian and normal to the 
shell meridian, may be seen in Figure 1 6 . These distortions are related 
by the equation 



6 = w cos <p - Ug sin <p (68) 

Therefore, the axial distortion at any point is directly related by the 

coordinates of the point and the deformation of the shell meridian. One 

finite difference solution for the linear behavior of shells of revolu- 

tion has been programed by gchaeffer for the digital computer and is 

presented in Reference 16. Using this program, a digital computer will 

calculate stress and moment resultants and displacements for thermal 

> •> ' 

and force loading varying along the meridian of the -shell* 

« r 

Distributed Thermal Inputs * 

1 * ^ ’I 

An analytical model of the telescope mirror was developed using 

f 

the finite difference solution of the shells equations from Reference 16. 

. * 

Utilizing this’ analysis, the distortions due to the application of a ‘ 
unit thermal input at the control point of interest was determined. 



kk 
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Assume the control point of interest is number 1. If we assign an 

amplitude A = 1 and TgjT^, . . . ,T n an amplitude A = 0, then 



By applying a unit amplitude at each control point separately, each 
column of coefficients may be determined. 

The analytical model used in this analysis has the capability of 
accommodating 502 control points along the shell meridian. A large 
number of control points is desirable because it increases the accuracy 
of the active optics system. However, since the flexibility matrix 
must be inverted in order to determine the stiffness matrix, the number 
of stations may be limited by the inversion routine. The inversion of 
large matrices is time consuming even for the best digital computers 
available. The accuracy of the inversion is limited by the behavior 
of the original matrix. One way to avoid this is to allow the thermal 
input to span several points in the finite difference analysis. The 
midpoint of the thermal input may be selected as the _ control point for 
correction by the figure error sensor in the active control system. 

The thermal inputs necessary for the correction of any given distortion 
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is therefore the amplitude (A) of the same distribution, used to deter- 
mine the flexibility matrix. 

Structures subjected to concentrated force loading have symmetri- 
cal flexibility matrices. Linear transformations with real symmetric 
matrices dominate the study of deformations of elastic media. Therefore, 
it is of interest to examine the flexibility matrix formed by distributed 
thermal inputs. Matrices generated by concentrated force loads are 
symmetrical because of the reciprocity theorem. This theorem states 
that for linear structures a force F^ acting through a displacement 
caused by force F-: does the same amount of work as force F- acting 

•through a displacement caused by force F^. This can he expressed in 
equation form as 


p i<°uV “ Wi) 


(70) 


and, therefore. 


c id = c ji fa> 

Since the coefficients of the flexibility matrix form a symmetrical 
array, the inverse or stiffness matrix must also he symmetrical. The 
basis of the reciprocity theorem lies in the fact that the total energy 
stored in an elastic system is independent of the order in which the 

loads are applied. This is also true of elastic systems deformed by 

/ 

thermal inputs. Therefore, the systems are analogous and for concen- 
trated thermal loads the flexibility and stiffness matrices are 
symmetric. 
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For distributed thermal inputs, however, the flexibility matrix 
may not be symmetric. In order to exa m i n e the flexibility matrix formed 
by distributed loads, we will now examine the case for a simply supported 
beam. A continuous beam with concentrated loads applied along three 
control points is shown in Figure 17(a). The deformation at any station 
along the beam may be determined from the following equations: 



If only the deformations under the loads are of interest, the equation 
can be simplified to (since T-^TyT^, etc. = 0) 




; 25 c 28 T 2 


c 55 c 58 ( t 5 

c 85 c 88 t 8 


(73) 


The application of a distributed load of unit amplitude on the 


same beam wide enough to cover several stations is shown in Figure 17(h). 
The deformations at control points 2, 5, and 8 due to the distributed 
load at stations 1, 2, and 3 are 




a. Beam with concentrated loads 


b. Beam with distributed loads 


£ 


Figure 17.- Sketch illustrating elastic beam with applied loading. 



49 


&2 - c 2]_(K) *** c 22^^ *** c 25( K ) 
6 S = c 51 (K) + c^p(l) + e^(K) 


52 v 


-55 \ 


5 

& 8 = c 8 l( K ) + c 82^4) + £35 (K) 


(74) 


j 


where K is a constant less than 1. If a distributed load of this 
form is always applied* the amplitude of the distribution may be denoted 
as T. If the deformations at points 2* 5^ und 8 are the only control 
points of interest* a new equation may be written as follows: 


f -1 
& 2 

f 

E C 2 l( K ) * c 22 ^ + c 23 ^ K H {_ 

— 

3 C 3 

( \ 
T 2 

> & 5 

> - 

[^051 (k) 4- 052 ( 1 ) + c 53 (k)J C 

L_J 
( ! 

1 / 

' T 

5 


[ 

[c 8 i(K) + cg 2 (l) -4 c 85 (K)] [ 

3 [ 3 

s 


The second column of coefficients for control point 5 due to the same 
load applied at stations 4* 5 -- and 6 is shown in the following equation: 


C x 
8 2 


jc.pi (K) + Cpp(l) 4 - 025 (K)j [p2i].(K) + 025(1) + C 25 (K)J £ 

(-■ 

T 2 

< 6 5 


[e 51 (K)+c 52 (l) + c 53 (K)3 [c^(K) 4 c 55 (l)+c 56 (K)] [ J 

1 5 

s 8 


(c 8 i(K) +c 82 ( 1 ) + <235 (K)j jjcg^CK) + c 85( 1 ) + c 86^ k 0 C 3 

s 


(75) 

The above equation indicates that the matrix of coefficients for the 
distributed load is not symmetrical. For example* 


= 51 (K) + C 52 (l) + c 53 (K) f c 2 i,(K) + e 25 (l) + e 2 g(K) 


(76) 
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because even though 


c 52 “ c 25 


( 77 ) 


the remaining coefficients 


C 51 + c 55 ? c 2 k + c 2 6 


( 78 ) 


are not necessarily equal. 



CHAPTER V 


RESULTS AND DISCUSSION 

The thermal Influence coefficients for the primary telescope 
mirror were determined using the computerized finite difference solution 
of the elastic shells equations. Ten control points equally spaced 
along the meridian of the primary telescope mirror were selected. The 
thermal inputs were considered to he applied by sources with controlled 
temperature distributions located on the back surface of the mirror. 

The thermal gradients in the interior of the mirror were neglected 
since they were shown previously to he quite small. 

The axial distortion due to the application of each symmetrical 
thermal input is shown in figure 18. These curves illustrate the dis~ 
tortion due to the application of a thermal input of unit amplitude at 
each control point. The location of ea,ch control point is indicated 
on the figure. The distortions are a maximum at the inner unrestrained 
boundary (except for control point 7) and decrease uniformly to zero 
at the axially restrained outer boundary. Maximum distortion is 
obtained by thermal inputs located at control points 4 and 10. It was 
noted previously that mirror distortions must be maintained to within 
2 microinches in order for the telescope to operate at the diffraction 
limit. The distortion created by a unit thermal input at control 
point 10 is over 300 times the allowable level. The large change in 
distortion between control points 9 and 10 as opposed to the relatively 
small change for points 3 and’ k indicates that this arrangement of 
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Axial distortion x 1(F (inches) 






Axial Distortion x 10^ (inches) 







control points may not be an optimum selection. For a thermal active 
optics control system, more control points at the outer region of the 
mirror should be investigated. Therefore, an examination of a control 
system based on equal annular areas as opposed to equal meridional 
increments may be desirable. 

The deflection influence coefficients were determined from the 
digital data used to plot the curves shown in Figure l 8 . By examining 
the distortion at all control points due bo a unit thermal input at each 
control point, each column of coefficients in the flexibility matrix was 
determined. These coefficients are shown in Table III. This matrix was 
inverted using the Jordan method which has been programed for the 
digital computer. The Jordan method is a library subroutine known as 
MATJ2TV and is available in the Langley program library. The program 
written to utilize this subroutine is listed in Appendix 3- The stiff- 
ness matrix resulting from this inversion of the flexibility matrix is 


given in Table IV. These coefficients give the amplitude of the thermal 
input necessary to correct a given set of axial distortions at the 
control points. In addition to the amplitude of the thermal input, we 
can also determine the accuracy of the amplitude of the thermal input 
necessary for diffraction-limited operation. This accuracy is given by 
the product of the minimum contribution (minimum coefficients K- - ) and 
the maximum tolerable error. An examination of the stiffness matrix 
indicates that the minimum contribution is 9 .I 9659 X 10 . The maximu m 


tolerable error for diffraction-limited operation was noted previously 
to be 2 microinches. Therefore, the amplitude of the thermal input must 
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lie controlled to less than 

A = 9-19659 X 10 +2 • 2 X 10“ 6 « 1.8 X lCf 5 °R 

In order to illustrate the use of the thermal active optics 
system, the distortions due to an acceleration-type loading were 
examined. The axial distortion due to a O.Olg static acceleration load 
was determined using the computer program of Reference 1 6 . The axial * 
distortion as a function of radius for this loading is shown in 
Figure 19- Also indicated on the figure is a sketch indicating the 
shell loading and positive direction of the distortion (5). For this 
relatively light loading, the axial distortion exceeds the tolerance 
limit indicated on the figure by a factor of about 25 . In order to 
correct this distortion, it will be necessary to introduce a distortion 
of equal magnitude and opposite sign by use of thermal inputs. The 
amplitude of the corrective thermal distortions were determined from 
the stiffness matrix of Table IV. These amplitudes were rounded to the 
nearest 0.001° R. To check the thermal inputs, the axial distortion 
due to both the acceleration load and the corrective thermal inputs 
were determined using the computer analysis of Reference 1 6 . The 
combined distortion is shown in Figure 20. Also indicated on the figure 
is the amplitude of the thermal input for each control point. The 
maximum residual distortion is well within the tolerance limit of 
2 microinches. It is significant that only very low amplitude inputs 
are necessary to correct the distortions which exceed the tolerance 








Control Point 


Amplitude of 
thermal input 


X 0.152 

2 0.193 

3 0.166 

4 0.171 

5 0.166 

6 0.177 

7 0.155 



Radius (Inches) 


Vn 

-1 


Figure 20.- Axial distortion of beryllium telescope mirror subjected 
to O.Olg acceleration type load and corrective thermal inputs. 
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limit by such a large amount. The feasibility of a thermal active 
optics system has thus been demonstrated for a cosine-type thermal 
input. 



CHAPTER VI 


CONCLUSIONS 

A thermal active optics technique to correct distortions in a 
thin telescope mirror has been developed. This technique utilizes 
measured distortions to determine the amplitude and location of thermal 
inputs necessary for correction of surface errors. 

The use of this technique has been demonstrated analytically by 
using a beryllium paraboloid. This mirror had 10 equally spaced con- 
trol points actuated by symmetrical thermal inputs. The stiffness 
matrix indicates that for this particular configuration the amplitude 
of the thermal inputs must be controlled to within 1.8 x 10”^ degrees. 
Distortions due to an acceleration-type loading were examined. These 
distortions exceeded the allowable tolerance by a factor of about 25 . 
Even though these distortions were quite large, they were corrected 
to well within the required accuracy. Axial temperature gradients 
resulting from the application of the thermal inputs to the rear 
surface of the mirror were shown to be quite small. 
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TABLE I. - PROPERTIES OF BERYLLIUM 


Property 

Value 

Modulus of elasticity (E) 

40 x 10 6 psi 

Poisson’s ratio (v) ' 

0.08 

Coefficient of thermal expansion (a) 

6.9 x 10“ 6 °R -1 

Thermal conductivity 

92 Btu/ft-hr-°R 

Density 

0.066 Ib/in^ 




I ABIE IE- CALCULATED VALUES OF THE VIEW FACTOR AND 


HEAT FLUX AT EACH CONTROL POINT 


Control 

point 

Radius * P™ 
(in.) 

- 

View 

factor* 

-mp 

Heat flux* Q/A 

(Btu/hr,in^) 

1 

6 

0.03366 

1.4229 X 10"? 

2 

11.4 

0.03384 

1.4243 x 10"3 

3 

16.8 

0.03400 

1.4250 X 10" 5 

1}- 

22.2 

0.03413 

1.4257 x io-3 

5 

27.6 

0.03424 

1,4264 x 10 “5 

6 

33-0 

0.03433 

1.4264 x io"3 

7 

38.4 

0,03440 

1.4271 x io“3 

8 

43.8 

0.03444 

1.4271 x io“3 

9 

49.2 

0.03447 

1.4278 x io"3 

10 

54.0 

0.03448 

1.4278 x io"3 




TABLE III 


FLEXIBILITY MATRIX GO 


■^TBIOas'H-SO 2.7U3aSE®01 3»2649e0E*01 3»343500C*01 

J ~ ' * ' " 

^^SOSOpaE^aa 244576iJ5E*01 3.951900E+0I 3.193600E*01 

Js;<i73aoaE*ag . 2aS2£a5£*01.2 B 7768OOE*Ql 2 = <?975Q0£>QI 

u 

|s*394gO0E«OO i«8H8SO0E®0l 2.4251O0E»01 2„723200E*01 
x jfcifc04OOBE»O0 1.5S36a8£»al 2.043000£*01 2.360600E»01 
' foiSPBPoaE« g a i.2lT9D0F«P 3 1. 66313 QE^OlI. U956S.0.QE*i!l 
^iol008E-»P0 9.445005E'*OO 1»29920DE*01 1,545900E«01 
L‘.96£9EDE-*00 6.909000E*0D 9.545O00E*OO l«i43100E*Oi 


E XB«>20E*0D 4.5260DOE*Q0 &.266000E*00 7.52B0OBE-»OO 
14CD00E- 03 2 .239 aOOE«OQ 3.102_000|>00 3,731 00QE«OO 


2c 9564Q0E»9 1 2jlO5l00E>01 7.76100gE*00-l .058000E*0 1-3.43Q500E+0 1-6.365200E+01 


2.B767G0E+G1 2 o 095S0GE*0I B.342000E*00-9.335000E*00-3.23B400E+01-6. 100900E*01 
2»768_60Q£»01 2*0759POE*01 9,OOBOOOE*00-7.821000E*00-3.000000£*0 1-5. 7691 OOE*01 
2.607200E«01 2,029000E*01 9.672000E*00-6.025000E*00-2.703000E+01-5.345200E*01 

r 


2.361900E*01 1.92990QE*01 
2.0195tt0£*01 a„744&QQ£*ai 
1,625600E*01 1»4567QOE»01 
1.214e00E*01. i;U0900E*01 
8.041Q0GE+00 7.42S0O0E*0O 
3»993000E*00_ 3,7<>Q0Q0E*00 


U014500E+01-4.059000E«00-2.349300E»01-4.619SOOE + 01 
l»O159OOE*Ql-2.il5OOQE»0Q“l .946200E*01-4«190400E*01 . 
9.354000E*00-4.590000E-01-I.S16400E*01-3.46130SE*01 
7.548000E*00 5.4IOOOOE-Ol-l.079200E*9l-2.6456oaE*Ol 
5.181000E*00 6.89Q000E-01-6. 744OOOE*0O-1 .76BB00E *01 
2.606000E*00 4. 050000£-01-3.260000E*00-8. 7470006*00 



TABLE 17.- STIFFNESS MATRIX jjK_ 


X 


1 5ilftZM;ME^V~l«2a640TE*0l 9.52l99ie*<|0-3.59Q353E*00 1.6076286*00-S.519662E-01 1.8OJ175E-01 8.6SF122E-Oi-3.7ttS10E*00 5.930679E*00 
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’ iJO 99776-02.-5. 363172E--01 2.6811696*00-7. 5920886*00 1.1930606*01-1.0566726*01 S.76Q626£*00-» .8975966*00-2.6673806*00 5.5319116*00 

1 F.3S0 1266-0 2 2.0676056- 02*8.8 867706-0 1 2.5890 686*0 0-7.5885886*00 1.1883666*01-1. 0962966*01 6.2618176*00-6. 2673616*00 7.0S3S90E*00 

1 _ __ ^ _ _ 

l 2U30 00 066 -0 2 - 7 *0 0 a 0 9 76 -0 2 J. 0981276-0 t r 8.2S_37896-01_ 2.8876386*00 7 7.S23|826»00 1. 181 7Sl£*0l -9.89852 76*00 1.97/6026*00 3.6887326*00 
' 2.9698106-02- 5.5110886-0 2 ' 1 .6631 36E-02 6 . 3522886-02- 3. 9598816-0 1 2,8188226*00-7.819858E*00 l. 2260716*01-1. 82822«E*0l 1.1388556*01 

l «... . 

" m08728t«02-2. 8551166-02 3. 2831216-02-3. 25S898£-02'V.SS79086-02-ii268788E-0 1 2.3283276*00-6.628832E*00 8.0221166*00-3.3710386*00 

~ 8 <02710 16 *02-6. 7895856-02 5.3838636-02 9.1965906-03-1.9281126-02-1.8929076-02-2.2502236-01 2.0696666*00-1.1291666*01 1.50881SE*01 | 

*— " ' ’ * „ “Tl 


not reproducible I 



APPENDIX A 


DIGITAL COMPUTER, PROGRAM TO GENERATE THE STEADY-STATE 


TEMPERATURE DISTRIBUTION IN A CIRCULAR DISC 


JOB. A030 1 » 1 .MARVIN RHODES ' » RD021 2*11 487201 1 

RUN ( S ) 

SETINDF. 

LGO. 

PROGRAM* MDR < lNPUT»0UTPUT»TAP£5rMNPUT,TAPE6=0arPUTf 

C THIS PROGRAM WILL GENERATE THE STEADY STATE TEMPERATURE 'DISTRIBUTION 
C IN A DISC OF THICKNESS T 

external func 

COMMON XLAM (20 ) » I 

DIMENSION D(1>»FOFX<1>»ANSJN(400>»ANSJ<400)»ANS<400)»F(20> 

C GENERATES 3(N> COEFFICIENTS DESIGNATED AS F(N) 

PRINT 3 

3 FORMAT ( » COEFF. NUMBER LAMBDA BESSELS FUNCTION CO£f 

IF. VALJE*) 

READ 1, A»B» N 

1 FORMAT (2FI0.5, 16) 

DO 7 1=1*20 

10 READ 2»XLAM<I) 

2 FORMAT (E16.8) 

CALL M GAUSS < A, B* N,D*FUNC»FOFX, 1 ) 

CALL BSSLS(XLAM<I)tANSJNtO,IERR) 

C=ANSJNi(l)»ANSJN<l> 

E=I.0*EXP<-XLAM(I)/25.) 

F < I ) =~2.*0 { 1 )/ (C*E) 

PRINT 4iI,XLAM<I)»ANSJNU)»B(I) 

4 FORMAT ( 5X»I3*3X*EI6.8,3X»E16.8*3X*E16.8) 

7 CONTINUE * • . 

C GENERATES TEMPERATURE AT 0.05 RADIAL INCREMENTS, FOR EAC H V ALUE O F ZETA 
30 IF (EOF,5)70»2Q 
20 READ 15, Z 
15 FORMAT (F10. 5) 

PRINT 17, Z 

17 FORMAT (3H Z=,F1 0 .5,2HXL) 

PRINT 18 ■ f 

18 FORMAT (6X*6HRADIUS»12X»HHTENPERATURE*/L 

R=0.0 ’ 4 

40 T=1 • 0 

DO 50 1=1,20 
G=XLAM<I)«R 

CALL BSSLS (G,ANS,0, IER) 


66 



H=EXP < XLAM < I ) «M Z-2 . ) /50 . ) ♦EXP < -XLAM < I > *Z/50 . 0 > 

HH=F(I)#ANS(1)*H 

t=t *hh 
50 CONTINUE 

PRINT 60»R*T 

60 FORMAT (616.8*5X»E16»8> 

R=R*0.05 

IF (R*LT # 1 • 05) GO TO AO 
GO TO 30 
70 STOP 
END 

SUBROUTINE FUNC(R,FOFX> 

DIMENSION FOFX(l)tANSJI400) 

COMMON XLAM (20 ) * I 
A=COS(6.283185307*R) 

b=xlam<ii»r 

CALL BSSLS(B,ANSJ*OtIERR> 

FOFX(l>=R*A*ANSJ(l> 

RETURN 

END 



APPENDIX B 


DIGITAL COMPUTER PROGRAM TO INVERT THE FLEXIBILITY MATRIX C 
TO OBTAIN THE STIFFNESS MATRIX K 


‘Job Vl» 0400% 450pbf. " A030 1T ' It MARVIN RHODES » RD02I2» H48»201I 

RUN (S) .1 "_J. 

SEIINDF. . . 

' P< PROGRAM. MAT- -aMPUT-»0UmU.T-.TAPE5=INRUT*TAP-E6=0UXRaT.)- 

-C--IHIS PROGRAM -MTLL—INVER-7— A~MATR‘IX— CALL-ED-A- BY SUBROUTINE— MAT-IN V. — iWPUl— DATA 
C ES -MR (.NUMBER OF. -R0MS3--.E0LL0ME0-~B.Y- -MATRIX.- A-OAXA—REAQ- B Y RO W S ^-INPU X- ANO— - 
OUTPUT FERMAT STATEMENTS (2..AND. 4 JL.MUSJ- BE. CHANGED .£0R„EACH. MATRIX.- - 

DIMENSION A ( 1 00 * 1001 t B_( 100*1 )_» IBI V 0I_U _Q O.L, .LND E XJ 100 ,2 _>. 

DIMENSION. DELTAjLIO L*.T (101 

READ 1» NR 

1 FORMAT (13) — - 

READ 2» KA(I*Dj!J=1-1lNR).».I=1jNR1™ „ . 

2 FORMAT ( (10F8.5L1 
PRINT 3 

3 FORMAT <//40X» ^ORIGINAL MAIRI**///}. 

PRINT 4*<<A(I,J)>J=1»NR)»I = 1»NR). - . 

CALL MAT1NV (A, NR »B* 0 .DETERM. IPIVOT • INDEX, 1 00 . ISCALE) 

I.._ PRINT 5‘ 

5 FORMAT(//540X*"INVERSE OF MATRIX A*///) 

PRINT ^t , HA(I»JltJ=l*NR) *1=1* NR) 

4 FORMAT ((10E13.6 , //)) 

DET=10.*»(100*ISCALE)*DETERM 
PRINT &» DET 

_6 FORMAT <///* VALUE OF DETERMINATE IS *E16.8> 

READ 7» (DELTA(I) »I=i»10> 

7 FORMATU0F8.6) 

DO 9 1=1(10 

T (11=0.0 
DO 8 J=l,10 

8 T<I)*A(I.J)*OELTA<J>.*T(I) 

PRINT10»I»T (I) 

10 FORMAT!* TEMP*»I2***«*E16.B> 

9 CONTINUE 
STOP 
END 
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